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Two-atom systems in small traps are of fundamental interest, first of all for understanding the 
role of interactions in degenerate cold gases and for the creation of quantum gates in quantum 
information processing with single-atom traps. One of the key quantities is the inelastic relaxation 
(decay) time when one of the atoms or both are in a higher hyperfine state. Here we measure 
this quantity in a heteronuclear system of s 'Rb and 85 Rb in a micro optical trap and demonstrate 
experimentally and theoretically the presence of both fast and slow relaxation processes, depending 
on the choice of the initial hyperfine states. The developed experimental method allows us to single 
out a particular relaxation process and, in this sense, our experiment is a ’’superclean platform” for 
collisional physics studies. Our results have also implications for engineering of quantum states via 
controlled collisions and creation of two-qubit quantum gates. 


Introduction 

The studies of two-atom systems in small traps at¬ 
tract a great deal of interest, in particular for engineer¬ 
ing of quantum states via controlled collisions and cre¬ 
ation of quantum gates in quantum information process¬ 
ing with a set of single-atom traps [1, 2]. The crucial 
points are the decoherence time and the lifetime related 
to the interaction-induced inelastic decay of a higher hy¬ 
perfine state. On the other hand, this type of inelastic 
processes, in particular heteronuclear ones, are important 
for the creation of multi-species quantum degenerate sys¬ 
tems [3], for obtaining ultracold heteronuclear molecules 
[4], and for ultracold chemistry [5]. However, the main 
obstacle in the studies of inelastic heteronuclear collisions 
in a trapped gas [6] is a simultaneous presence of a large 
variety of loss mechanisms, which complicates the analy¬ 
sis. In magneto-optical traps where many heteronuclear 
systems have been studied [7-12], aside from collisional 
processes one has radiative escape. In optical dipole traps 
there are homonuclear inelastic collisions [13-15], and at 
sufficiently large densities three-body recombination be¬ 
comes important [16]. It is therefore crucial to perform 
experiments allowing one to single out a particular in¬ 
elastic process [17-19]. 

This is done in the present paper. We study a two- 
atom system of different isotopes of rubidium (single 
85 Rb and single s 'Rb) in a micro optical trap. One of 
them or both are in a higher hyperfine state, and we 
measure the corresponding rate of inelastic relaxation ac¬ 
companied by the loss of the atoms. The homonuclear 
collisions are absent and our measurements give pure loss 


rates of specific hyperfine heteronuclear collisions. The 
experiments are done at temperatures close to the bor¬ 
der of the ultracold limit (tens of microkelvins) and are 
supported by finite temperature coupled channels calcu¬ 
lations. Our work can be easily extended to other alkali 
atoms, even to atom-molecule collisions [20, 21], thus al¬ 
lowing further understanding of heteronuclear collisions, 
a precise test of atomic collisional theory, and applica¬ 
tions to quantum information processing. 


Experimental setup and results 

Our two-atom heteronuclear system is composed of a 
single 85 Rb and a single 87 Rb in a micro optical dipole 
trap (ODT), and there are three important points in the 
experiment. The first one is a sequential trapping of 
a single 87 Rb in a static ODT and a single 85 Rb in a 
movable ODT [22], and we make sure that two atoms 
of different isotopes are actually trapped (see Fig.l and 
Methods). Second, we shift the movable ODT to over¬ 
lap with the static one, and adiabatically turn off the 
movable trap. We get 87 Rb and 85 Rb in one trap with 
probability of about 95%. The third point is that because 
of collisional blockade [23] we have to kick out one of the 
atoms before detecting the presence of the other one. By 
optimizing this procedure we have minimized unwanted 
atom losses to less than 3%. 

Depending on the hyperfine states of 87 Rb and 85 Rb, 
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FIG. 1. (color online) Experimental setup and measurement time sequence, (a) Schematic diagram of the experimental setup. 
Two 830 nm lasers are collimated and then strongly focused by an objective (Linos,HALO30) into the vacuum chamber to form 
two ODTs. The movable ODT is from 830 nm laser-1 and can be shifted to overlap with the static ODT (from 830 nm laser-2) 
by controlling PZT-1 (Piezoelectric Ceramic Transducer). The fluorescence of trapped single atoms is collected by the same 
objective and guided to SPCM (Single Photon Counting Module) for detection. PZT-2 controls the fluorescence collecting 
region. A detailed description can be found in Methods. (b)Time sequence in the experiment. Each survival probability in our 
experiment is the result from 300 repeated measurements. 


there are three inelastic decay processes: 

(A) 87 Rb(F = 2) + 85 Rb(F = 3) =► 

87 Rb(F = 1) + 85 Rb(F = 3) 

87 Rb(F = 2) + 85 Rb(F = 2) 

87 Rb(F = 1) + 85 Rb(F = 2) 

(B) 87 Rb(F = 2) + 85 Rb(F = 2) => 

87 Rb(F = 1) + 85 Rb(F = 3) 

87 Rb(F = 1) + 85 Rb(F = 2) 

(C) 87 Rb(F = 1) + 85 Rb(F = 3) => 

87 Rb(F = 1) + 85 Rb(F = 2) 

We have not set a magnetic field, and for each atomic 
spin F in the initial state of the collision the states with 
all possible values of the spin projection Mp are likely 
equally populated. The energy released in the inelastic 
processes A, B and C is about several GHz and it exceeds 
the trap depth Uq by more than two orders of magnitude. 
Therefore, both atoms are ejected from the trap as a 
result of the inelastic relaxation, which is confirmed in 
the experiment by selectively kicking out 85 Rb or 87 Rb in 
the B process. In most of our experiments the trap depth 
is Uq = 0.6 mK, which in our configuration provides the 
radial trap frequency = 38.8 ± 0.1 kHz and the axial 
trap frequency lo z = 3.2 ± 0.1 kHz. 

We measure the survival probability P(t) for the atoms 
to remain in the trap at time t (see Fig. 2). For each t we 
execute 300 repetitions of the loop sequence of Fig. 1(b). 
In the case of A and B processes the decay is strongly 
dominated by the interaction-induced spin relaxation. 
The probability P(t) is then described by an exponential 
time dependence. Within less than 10% of uncertainty 


the experimental data can be fitted with an exponential 
function P = wexp(—t/r) + wq. The presence of the 
residue wq has two reasons. First of all, the two-atom 
system is obtained with 95% probability, and there are 
traps with only one atom which remains trapped on a 
much longer time scale (about 11s [24]) than the col- 
lisional lifetime r. Second, doubly polarized pairs (for 
each atom the spin projection is equal to the spin) can 
decay only due to weak spin-spin or spin-orbit interac¬ 
tions and practically remain stable on the time scale of 
our experiment. 

The process C is much slower than A and B, and our 
measurements for this process have been made on a time 
scale of about 500 ms. In this case the decay is signifi¬ 
cantly influenced by single atom spin relaxation, and we 
have to take it into account in the rate equations for ex¬ 
tracting r from our measurements (see Methods). 

The measured r has about 15% of statistical uncer¬ 
tainty which decreases with increasing the executed loop 
numbers. Aside from single atom spin relaxation, the 
value of r is influenced by single atom loss events. The 
heating rate in the dipole trap is about 20 fiK/s [25] and 
it increases the collisional volume, thus slightly increas¬ 
ing the decay time r. We estimate the overall uncertainty 
in our values of r as about 20 %. 

For the A collisional process, we also vary the temper¬ 
ature in order to test the dependence of r on the effec¬ 
tive volume (density) of atoms in the trap. As expected, 
the time r increases with temperature and one can see 
this from the comparison of the results in Fig. 2(b) and 
Fig. 3(a). We then tested that the result for r does not de¬ 
pend on whether we kick out 85 Rb or 87 Rb for measuring 
P(t) (see Fig. 2(c) and Fig. 3(b)). It is easy to conclude 
that not only relaxation times are very close to each other 
(as one sees in Table 1), but also the functions P(t). 
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FIG. 2. (color online) Experimental data for the decay rates, (a) Energy levels of hyperfine states of 87 Rb and 85 Rb. (b), 
(c), and (d) Survival probability P versus time t for the A, B, and C collisions, respectively. The measurements are done for 
the survival probability of 8 'Rb after kicking out 85 Rb. The black squares are experimental data, with each point being the 
result from 300 repeated measurements. In b) and c) the solid curves show a fit by the formula P = wexp(-f/r) + wo, and in 
d) a fit with the numerical solution of the rate equations including single atom spin relaxation. The data are collected at the 
trap depth Uo = 0.6 mK, and the initial temperatures Ts 7 = 35 ± 3 pK, Tss = 15 ± 1 pK. Numbers in bracket are the fitting 
standard deviations. 



FIG. 3. (color online) The decay under different conditions, (a) The same as in Fig. 2(b) but for Ts 7 = 47±3 pK and Tss = 27±2 
fj, K. (b) The same as in Fig. 2(c) but for kicking out 87 Rb. 


Theory and analysis 


The rate equations for the inelastic decay processes A, 
B, and C can be written as 


dP _ P 
dt t 


( 1 ) 


where P(t) is the probability that at time t the atoms are 
still present in the trap, and r is the relaxation (decay) 
time that we measure. These processes occur at inter¬ 


atomic distances of the order of or smaller than the radius 
of the interaction potential R e = (mCe/H 2 ) 1 ^ 4 ~ 80A 
(Cq is the Van der Waals constant). At our trap fre¬ 
quencies and temperatures from 15 to 55 /xK we have 
T fko p ,huj z , and the motion of atoms in the trap is 
surely quasiclassical. The extention of the wavefunction 
of the atoms in the trap is rx ~ ( T/mui 2 ) 1 / 2 ~ 0.2/im, 
where m and u> are characteristic values of the atom mass 
and trap frequency. Thus, we have the inequality 

r T > R e - 


( 2 ) 
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Therefore, the decay (relaxation) time r can be expressed 
through the relaxation rate constant a in free space: 


1 a 

T = V&' 


(3) 


where 




/ \ 3 / 2 
/ 27rT e ff \ 

l 4/3 2/3 1 


(4) 


is the effective volume, T e g = p(Ti/mi + T 2 /m 2 ) is the 
effective temperature, with mi,m 2 and T\,T 2 being the 
masses and temperatures of 87 Rb and 85 Rb, and p is the 
reduced mass. The rate constant a in Eq.(3) is averaged 
over the Boltzmann distribution of relative momenta k 
at effective temperature T e g 


/ h 2 N 3/2.00 

V ) Jo 


( k )exp 


h 2 k 2 \ 

2/tTeff ) 


47t k 2 dk. 


(5) 


Note that due to elastic collisions between the atoms the 
two-atom system eventually acquires an equilibrium tem¬ 
perature (Ti +T 2 )/2. However, it is different from the 
initial effective temperature T e g by less than 1%, and so 
will be the effective volume and the average value of a. 

At our temperatures the quantity kxR e ~ 0.5, where 
kx = ( mT/h 2 ) 1 / 2 is the thermal momentum, so that we 
are close to the border of the s-wave scattering limit. 
Therefore, in addition to the s-wave scattering, we took 
into account the scattering with higher orbital angular 
momenta. The rate constants a a, &b, and ac for the 
processes A, B, and C were calculated using the cou¬ 
pled channels method [4, 26] at finite collision energies 
(see Methods). In the center of mass reference frame the 
Hamiltonian governing the collisions has the form 

H = ] k + ^ + Fel(r) + yhf ’ (6) 

where r is the interatomic distance, p is its conjugate 
momentum, and i is the orbital angular momentum op¬ 
erator. The interatomic interaction operator is given by 
V el (r) = V s (r)P s + V t (r)P t , with P s and P t being projec¬ 
tors onto the electronic singlet and triplet states of the 
colliding pair of atoms. The term Vhf = aiSi-Ii+a 2 S 2 -I 2 
is the hyperfine interaction, where Si, Ii, 01 and S 2 , 12 , a 2 
are the electron and nuclear spin operators and hyperfine 
constants for 87 Rb and 85 Rb, respectively. The total spin 
operators of the atoms are Fi = S 1 +I 1 and F 2 = S 2 +I 2 , 
and the total spin operator of the pair is F = Fi + F 2 . 
The Hamiltonian H of Eq.(6) conserves both the total 
spin F and its projection Mf- It also conserves the or¬ 
bital angular momentum t and its projection Me. 

The rate of inelastic spin relaxation occuring when at 
least one of the colliding atoms is in a higher hyperfine 
state, can be expressed through the real and imaginary 


parts of the elastic scattering amplitude. The well-known 
formula for the inelastic rate constant [27], averaged over 
the initial spin projections, is transformed to (see Meth¬ 
ods): 


a(F 1 ,F 2 ,k) = 


47rfi 


F—F\-\-F2 

B “+ 1 ) £ 


(2Fi + 1)(2F 2 + l)n _ 
x(2 F + l)[Im fe(F 1 ,F 2 , F, k) - k\fe(F u F 2 ,F, fc)| 2 ]- (7) 


Accordingly, a a = a(2,3), as = cx( 2,2), and ac = 
a(l, 3). The quantity fe(Fi,F 2 , F 7 k) in Eq.(7) is the am¬ 
plitude of elastic £-wave scattering of these atoms at the 
total spin F. 

We apply the accumulated- phase method (see Meth¬ 
ods) and calculate the accumulated-phase parameters 
from the known properties of homonuclear 87 Rb 87 Rb and 
85 Rb 85 Rb collisions using mass scaling [26]. The main in¬ 
accuracy of our calculations stems from the choice of the 
accumulated phase, and we have checked that our results 
are stable within 5% when the value of this phase is var¬ 
ied by a few percent. The p -wave (1=1) contribution at 
our temperatures is comparable with the s-wave (£ = 0) 
one, but the contributions of the d- wave and higher par¬ 
tial waves are below 1%. Therefore, in the following we 
confine ourselves only to the s-wave and p -wave scatter¬ 
ing. 

In Table 1 the results of the calculations for the pro¬ 
cesses A, B and C are compared with the experimental 
data. For the fast A and B processes one sees an agree¬ 
ment between experiment and theory within the error 
bars of the experimental data. For the slow C process 
the calculated a is near the upper bound of the experi¬ 
mental value accounting only for statistical uncertainties. 
The reason for this small discrepancy is that the heating 
effect, although fairly small, still increases the effective 
volume so that the measured r actually corresponds to 
slightly higher temperatures than the initial ones. This 
means that the experimental value of a at the initial 
temperatures should actually be slightly (by about 15%) 
higher than the one in Table 1. 


Discussion 

The relaxation rates obtained in our work are rather 
high. In interesting experiments with spinor heteronu- 
clear mixtures [28], this places an upper limit of at 
about n ~ 10 12 cm -3 on the density if a higher hy¬ 
perfine atomic state is involved. Our relaxation rate 
constants are several orders of magnitude higher than 
the ones for doubly polarized atomic states, like for ex¬ 
ample 87 Rb(Fi = 2, Mi = 2), where the relaxation is 
caused only by weak spin-spin and spin-orbit interac¬ 
tions [29]. We expect a similar reduction for inelastic 
collisions of doubly polarized 87 Rb and 85 Rb atoms, such 
as 87 Rb(Fi = 1,Mi = l)+ 85 Rb(F 2 =3, M 2 = 3). 
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TABLE I. Summary of the results (the experimental data are for the trap depth Uo = 0.6 mK (cj p = 38.8 kHz, to z = 3.2 kHz)) 


Collisional 

process 

Temperature and 
trap depth(pK ) 

Decay time (ms) 

Experimental a(cm 3 /s) 

Calculated a(cm 3 /s) 


T 87 = 35 ± 3, T 85 = 15 ± 1, kick out 85 Rb 

ta=74±12 

(5.9±l.l)x 10” 11 

5.6x KT 11 

A 

T 87 = 47 ± 3, T 85 = 27 ± 2, kick out 85 Rb 

ta=121±13 

(6.5±0.8)x HT 11 

5.9x KT 11 


T 87 = 35 ± 3, T 85 = 15 ± 1, kick out 85 Rb 

t b = 70±8 

(6.3±0.9)x HT 11 


B 

T 87 = 35 ± 3, T 85 = 15 ± 1, kick out 87 Rb 

ts=67±9 

(6.6±1.0)x HT 11 

6.8x KT 11 

C 

T 87 = 35 ± 3, T 85 = 15 ± 1, kick out 85 Rb 

tc =1800±250 

(2.4±0.4)x 1(T 12 

3.2x 10 -12 


Our results may also have implications for engineering 
of quantum states and creation of two-qubit quantum 
gates in a system of single-atom traps, in particular with 
respect to a proper selection of hyperfine atomic states. 
The relaxation time r in our experiment is still much 
larger than the characteristic time on which the single¬ 
atom coherence is preserved - the decoherence time Td c - 
The time Td c measured at the same temperature T is 
about 1 ms and it increases with decreasing T [30, 31]. 
On the other hand, the time r strongly decreases with 
T due to a related decrease in the effective volume V e g. 
The expected reduction in r for the atoms in the ground 
vibrational state is about a factor of 30 compared to the 
data in Table 1, and it becomes of the same order of 
magnitude as Td c , at least for the atoms in the states 
F\ = 2, F 2 = 2 and F\ = 2, F 2 = 3 (processes A and B). 
Even for the slowest process C the expected r is ~ 50ms, 
which is close to the maximum Td c obtained up to now 
in experiments with single-atom traps [32] . 

The time r can be certainly increased by decreasing the 
trap frequencies. However, there is a lower bound for the 
frequencies. This is because for quantum gates the time 
Td c and, hence, r should be at least 4 orders of magnitude 
larger than the operational time r op [33] . The latter is on 
the millisecond level in the considered situation and can 
not be much smaller than the inverse interaction energy 
of the two atoms in the trap, so that r op ~ hV e s/g, where 
g = 4irh 2 a/m is the coupling constant for the elastic 
interaction, and a is the scattering length. We thus have 
the condition: 

r > Tdc > 10 4 r op ~ (8) 

g 

As the effective volume is V e g oc 1/w 3 , for realistic T dc 
below Is the lower bound for the trap frequency is about 
a few kiloHertz. 

There is also a fundamental limitation independent of 
the trap frequencies. Since r = V e g/a, equation (8) im¬ 
mediately leads to the inequality 

■§- » 10 4 - (9) 

ha 

For common values of g/K from 10 -9 to 10 -11 cm 3 /s 
equation (9) can be satisfied for doubly polarized atomic 


states, such as 87 Rb(E = 2 ,Mp = 2) where g/h ~ 10“ 10 
cm 3 /s, and a < 10 -14 cm 3 /s. Another option would be 
to increase g by using a Feshbach resonance, although 
this can also increase inelastic losses. 

The advantage of our work is that we study collisions 
in a two-atom system, which allows us to single out a 
particular collisional process. In upcoming experiments 
we intend to prepare atomic states with given spin pro¬ 
jections Mf and cool the atoms down to the ground vi¬ 
brational state [34, 35] in order to execute a possibility of 
creating a quantum gate. Our work also paves a way to 
the creation of single heteronuclear molecules and to the 
studies of atom-molecule and molecule-molecule binary 
systems. 


Methods 

Details of experimental methods 
In our experiment single atoms are trapped in micro op¬ 
tical dipole traps (ODT) which are formed by strongly 
focusing 830 nm lasers to a beam waist of about 2.1 /im. 
The dipole lasers follow the paths shown in Fig. 1(a) by 
the red solid lines. The movable ODT is initially 5 gm. 
away, and it can be shifted to overlap with the static ODT 
by changing the voltage of PZT-1. We detect trapped 
atoms by collecting the fluorescence in the trap region as 
shown by blue dashed lines. The fluorescence is coupled 
into a polarization-maintaining (PM) single mode fiber 
(with the core diameter of 5 /irn) for spatial filtering and 
is then guided to a single photon counting mode (SPCM, 
AQRH-14-FC). Owing to a moderate core diameter of 
the fiber and to a fairly large distance between the two 
traps, we can selectively collect the fluorescence from one 
of them. The disturbance from the other trap (crosstalk 
effect) is eliminated by properly adjusting the voltage of 
PZT-2. 

The whole experiment is executed following the time 
sequence shown in Fig. 1(b). We first trap 87 Rb in the 
static ODT and then 85 Rb in the movable trap. We do an 
extra check of the 8 ' Rb presence in the trap by turning on 
the 87 Rb MOT again for 20ms at step 3. Only when 87 Rb 
is still trapped, we record the final result of the collision. 
In order to check that only a 8 'Rb is in the static trap 
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and meanwhile only a 85 Rb in the movable trap, we first 
detect the absence of the fluorescence of 8 ' Rb from the 
movable trap and the absence of fluorescence of 85 Rb 
from the static trap. This is followed by the fluorescence 
detection of 87 Rb and 85 Rb in the static and movable 
traps, respectively. 

In step 4 we first prepare 85 Rb in F 2 = 2 and s 'Rb 
in Fi = 1 states in order to eliminate the unwanted col- 
lisional losses when switching off the MOT repumping 
lasers 1ms before the MOT cooling lasers. For minimiz¬ 
ing the heating effect we have optimized the process of 
transferring 85 Rb to the static trap. The transfer effi¬ 
ciency can go up to 98% and is limited by the detection 
efficiency and heating losses. The probability that 8 'Rb 
survives when the movable trap approaches the static one 
and then is adiabatically switched off, is also about 98%. 
The temperature is measured in this type of process by 
using the release and recapture technique [24] when one 
of the traps (either static or movable) is empty. 

In this way we create a heteronuclear two-atom sys¬ 
tem. After a certain time, we kick out one of the atoms 
from the trap by using resonant lasers. Optimizing the 
laser intensities and shortening the pulse duration to 
0.1ms, we have minimized the unwanted losses to less 
than 3%. Eventually, we have succeeded in trapping 
two heteronuclear single atoms in the static ODT with 
about 95% probability. 

Analysis of the experimental data 
For the slow process C we extracted r from the fit of 
our experimental data with the numerical solution of the 
rate equations taking into account single atom spin re¬ 
laxation. These are linear equations for the quantities 
Pfi F 2 (t) representing the probabilities that both 87 Rb 
with spin i 7 ) and 85 Rb with spin F 2 are present in the 
trap at time t: 

dPi3 _ -Pi 3 (-P23 + -P12 ~ 2P13) 

dt TC T r 

dP\2 _ (-P13 + P22 ~ 2 Pi 2 ) 

dt T r 

dP23 _ P23 {P13 + P22 ~ 2F23) 

dt TA T r 

dP 2 2 _ P22 (-P23 + -P12 — 2P22) 

dt TB T r 

For the time of single atom spin relaxation we did an 
independent measurement, and the measured value is 
T r = 1100 ± 150 ms. The interaction-induced relaxation 
times ta and tb were taken from the data for the A and 
B processes. In the main text, Fig. 2(d) shows that the 
best numerical fit Pi 3 {t) deviates from the experimental 
data by less than 10%. 

Calculation of the inelastic rate constants 
At low temperature the inelastic spin relaxation occurs 


when at least one of the colliding atoms is in a higher 
hyperhne state. The rate constant of this process is given 
by [27] 


a(F 1: F 2 ,k) 


7 Th v y. (2£ + l)|5g ) | 2 

ak A —1 (2 Fi + 1)(2F 2 + 1) ’ 


( 10 ) 


The quantity S^' is the S'-matrix element for the Owave 
scattering from the initial state i characterized by the 
atom spins F\, F 2 and their projections Mi, M 2 to a final 
state / which has a lower internal (hyperfine) energy, so 
that there is an energy release in the inelastic scattering 
process. In Eq.(10) we also averaged over the initial spin 
projections Mi and M 2 . Due to the unitarity condition 
for the P-matrix elements we have: 

£ | S <«| 2 = 1 H *®! 2 - ( 11 ) 

f i'Ai 


where Su is the S'-matrix element for elastic scattering 
in which the spin projections Mi,M 2 remain the same, 
and Su> is the S'-matrix element for elastic scattering 
which changes Mi,M 2 to Expressions for the 

S-matrix elements through the corresponding scattering 
amplitudes read [27]: 

S$ =Sw+2ikf$(k). (12) 


The amplitudes fa/ are conveniently expressed through 
the amplitudes fe(Fi, F 2 , F, k) of elastic scattering of 
atoms with spins Fi,F 2 at the total spin F of the pair: 

F=Fi~\~F2 

f${k) = J2 (Pi Mi F 2 M 2 \FiF 2 FM) 

F=|Fi-F 2 | 

(F 1 M' 1 F 2 M' 2 \F l F2FM)h{F l ,F2, F, k). (13) 

In the absence of a magnetic field the amplitudes 
fefFi, F 2 , F, k) do not depend on the total spin projec¬ 
tion M. The Clebsch-Gordan coefficients which appear 
in Eq.(13) satisfy the summation rule: 

Y (FiMiF 2 M 2 \F\F 2 FM)(F 1 M 1 F 2 M 2 \F 1 F 2 FM) 

Mi,M 2 

= ^ff^mm- (14) 

Substituting Eqs.(12) and (13) into Eq.(10), and making 
use of Eq.(14), we arrive at equation (7) of the main text: 


a(F 1 ,F 2 ,k) 


1 Airh 

(2Fi + 1)(2F 2 + 1) IT 


E( 2£+1 ) 


^=0 


F 1 +F 2 

Y (2F+1) [Im f e (Fi, F 2 , F, k)-k\ f t (Fi, F 2 , F, k)\ 2 \. 

F=|Fi —F 2 | 


The amplitudes fe(F±, F 2 , F, k) were calculated numer¬ 
ically using the coupled channels method [4] . Our imple¬ 
mentation of this method is described in Ref. [36]. The 
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asymptotic behavior of the scattering states is enforced 
at a distance of lOOOao, with ag being the Bohr radius. 
The accumulated-phase boundary condition is applied at 
r 0 = 16a 0 . It summarizes the short-range physics in 
the region of distances r < ro, where the triplet and sin¬ 
glet interaction potentials are poorly known, into 6 phase 
parameters. We calculate these parameters for heteronu- 
clear 87 Rb 85 Rb collisions starting from the known data 
for homonuclear 87 Rb 87 Rb and 85 Rb 85 Rb collisions and 
using the mass scaling technique [26], which exploits the 
fact that the Born-Oppenheimer electronic potentials V s 
and Vt do not depend on the type of isotopes. The hy- 
perfine coefficients a\ and a 2 were taken from Ref. [37]. 
Solving the coupled differential equations for the wave- 
functions within the subspaces characterized by the con¬ 
served quantum numbers F, M we calculated all scatter¬ 
ing amplitudes .F^-F. k) as functions of the inci¬ 
dent relative momentum k. The calculated rate constants 
a(Fi , F' 2 , k ) are then averaged over the thermal distribu¬ 
tion of k according to Eq.(5). Note that our zero tem¬ 
perature result ac = 0.8 x 10“ 12 cm 3 /s is fairly close to 
the lower bound of the interval (1.2 — 4.5) x 10~ 12 cm 3 /s 
found for the C process with 87 Rb(.Fj = l,Mi = —1) 
and 85 Rb(F 2 = 3,M 2 = 3) in Ref. [38] which used old 
data for the interaction potentials. 
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